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Abstract 

We consider the motion of red blood cells and other non-spherical microcapsules dilutely sus- 
pended in a simple shear flow. Our analysis indicates that depending on the viscosity, membrane 
elasticity, geometry and shear rate, the particle exhibits either tumbling, tank-treading of the mem- 
brane about the viscous interior with periodic oscillations of the orientation angle, or intermittent 
behavior in which the two modes occur alternately. For red blood cells, we compute the complete 
phase diagram and identify a novel tank-treading-to-tumbling transition at low shear rates. Ob- 
servations of such motions coupled with our theoretical framework may provide a sensitive means 
of assessing capsule properties. 
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Human red blood cells suspended at low concentrations in steady simple shear flow have 
been observed to exhibit two types of motion. While cells suspended in a low viscosity 
media tumble continuously, cells suspended in a fluid with sufficiently high viscosity exhibit 
'tank-treading' motion P]. Here, we use the term tank-treading to describe a cell that main- 
tains almost constant shape and orientation in the laboratory frame, but whose membrane 
circulates around the interior much like the motion of a tank tread. Keller and Skalak 
analyzed the motion of a tank-treading viscous ellipsoid and concluded that the behavior 
depends on the ratio of the viscosities of the inner and outer fluids and was independent of 
the shear rate. However, it has recently been observed that red blood cells oscillate about a 
: ixed angle while tank-treading and that tumbling can be induced by lowering the shear rate 
^. Furthermore, the dynamics of synthetic microcapsules in simple shear flow have been 
shown to depend not only on the viscosity ratio, but also on the shear rate: at high shear 
rates, the capsule surface tank-treads about the interior while the orientation oscillates; at 
low shear rates, the capsule tumbles . These observations are unaccounted for by cur- 
rent theory and motivate the present work. In this letter, we provide a unified theoretical 
framework for analyzing the motion of both non-spherical capsules and red blood cells in 
simple shear flow and attribute the observed shear rate dependent behavior to a periodic 
variation in the elastic energy during tank-treading. 

Here, we model the cell or capsule as an ellipsoid of viscosity fi' contained in an elastic 
membrane and immersed in a fluid of viscosity fi. In the frame of the ellipsoid having axes 
of length Oj, the surface is at the position ^ = 1. The external flow unperturbed by the 
ellipsoid is a simple shear flow of rate 7 in the laboratory reference frame (see figure la). 
We assume that material elements in the membrane move along elliptical paths in planes 
parallel to the X1-X2 plane, with velocity field as in Then, the position of each membrane 
element is defined by a phase angle (j){t) and the tank-treading frequency is dt4>{t). When 
= 27r a point on the surface has returned to its starting point. Our formulation extends 
the analysis of for a viscous ellipsoid to include the effect of an elastic membrane and we 
refer the reader to the literature for the derivations of several equations unrelated to elastic 
membrane. Torque balance yields an equation for the rate of change in the orientation of 
the ellipsoid, 6, in the laboratory frame [2] 



dtO = --- ^—^dt(p+ -^—^ cos 29 1) 
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The tank-treading frequency, 9^0, is found by balancing the rate of work done by the external 
fluid on the ellipsoid with the internal dissipation and the rate of change of the elastic energy 
stored in the membrane. The rate of work done by the external fluid on the ellipsoid, Wp, 
and the internal dissipation, D, are P] 

Wp = Vii[f2{dt(j)f + /379i0cos2^], (2) 
D = V/i7l(5^0)^ (3) 

where V is the volume of the ellipsoid and fi are functions of the ellipsoid geometry defined 
in the endnote 

In a series of recent experiments on red blood cells in a transient shear flow, Fischer 
showed that after relaxation from deformation the rim of a red blood cell is always formed by 
the same part of the membrane. This shape memory could not be eliminated by continuous 
deformation for periods up to 4 hours, showing that the effect is not due to visco-elastic 
relaxation. These results imply the existence of an elastic energy that has a minimum when 
the membrane is in static equilibrium. Since an element displaced from the rim can return to 
either side of the rim depending upon its proximity 6], the elastic energy must pass through 
two minima during a tank-treading revolution. Similarly, for the case of a polyamide capsule, 
the unperturbed shape is a slightly non-spherical ellipsoid and we assume that the elastic 
energy during tank-treading has the same vr-rotational symmetry as the undeformed shape. 
Consequently, we take elastic energy to have the form 

E = Eosm\(f)). (4) 

Fischer's experiments IgJ can also be used to estimate the size of the elastic energy change, 
Eq. When the shear flow is stopped, the rate of work done by the external fluid Wp vanishes. 
The slender geometry of the red cell allows us to neglect dissipation in the external fluid 
in this analysis so that the rate of change in the elastic energy balances the internal 
dissipation (jH)) so that 

= Eosm{2(P)dt(j) + Vfi'h{dt<f)y. 

The volume of a red blood cell is V ~ 7 ■ 10~^''m^ and /i ~ 15/4. The internal viscosity 
fi' is due to a combination of the cytoplasm viscosity and the membrane viscosity. The 



dissipation due to the membrane viscosity has 



3een shown experimentally to be between 2 



and 4 times the dissipation in the cj^oplasm 



0,3, 



and we choose the value fi' = 4:fj,interior 



4- 10~^Pa-s corresponding a membrane dissipation equal to three times the dissipation in the 
cytoplasm. Fischer's data show the maximum tank-treading frequency to be approximately 
dt(j) ~ 0.25s"^ so that 

^0 = Wil-^t^Uax ~10-^'J. (5) 

We can now estimate the amount of strain corresponding to the stored energy Eq. The 
properties of the red blood cell membrane are determined by both the phospholipid bilayer 
and the adjacent protein network. The high resistance to change in membrane surface area 
insures that shear (area preserving) deformations dominate. The resistance to shear and the 

n 

shape memory are most likely due to the underlying protein network |6||. For a membrane 
of shear modulus E^^Q ■ 10~^N/m covering a surface 5* ~ 1 ■ 10~^''m^ and undergoing 
a shear strain u the elastic energy is approximately Eq ^ ^"'^^ . Equation Q then yields 
the characteristic shear strain u ~ 0.2, a reasonable value. 

Equating the rate of work by the external fluid (0) with the internal dissipation (jSj) and 
the rate of elastic energy storage (jH) yields 

Vij{f2dt(j)^ + Mdt(f) cos2e) = V^x'h{dt<j)f + Eo^m{2(t>)dt<i>, (6) 

which can be solved for the tank-treading frequency 

dt<P = T-^^ArriUe sin 20 - cos 29) , (7) 

where A = /i'/ /i is the viscosity ratio and f/g = is the ratio of the change in the elastic 

energy to the work done by the external fluid during a rotation, i.e., how stiff the cell or 
capsule is relative to the forcing from the external shear flow. Scaling time, t = T/j, yields 
the dimensionless system of equations for 6 and 

= r ^\ r (Ue sin 20 - COS 29) , 

Ore = -'-- 4^a^0+ i4f4eos2^. (8) 
2 ai + a2 2ai + a2 

The system's behavior is now completely determined by Ue-, A, and the geometric ratios 

ai/a2, 01/03. Although observations of microcapsules have shown that the capsule geometry 



periodic changes in axes length have little 



oscillates periodically during tank-treading 
impact on solutions of 

A simple way to understand the relationship between the tank-treading and oscillation 
periods is to expand (jH)) using Ue as a small parameter: = 0o + t/e0i5 6* = 6*0 + f^e^i- 



The 0{U^) solution recovers the case studied by Keller and Skalak where the tumbling 
transition is induced via an increase in the viscosity ratio and the steady solution is described 
by the constant orientation angle 9q and the tank-treading frequency (9^00 = ^- The 0{Ul) 
perturbation to this steady state is 9i = csin(2ti;T — 6), where the constant c and the phase 
lag 6 between the elastic energy and the orientation angle depend on the specific parameters 
of the system. Hence, 6 = 6o + UeCsm{2ujT — 6) and the two minima of the elastic energy in 
each tank-treading cycle cause the system to oscillate twice about a fixed angle in the time it 
takes a surface element to tank-tread back to its initial position. This effect is illustrated in 
figures lb and Ic and solutions of dHl) show that it is also true for large amplitude oscillations. 

nn 

Slightly non-spherical polyamide microcapsules |5[ and numerical simulations of slightly 
non-spherical liquid capsules enclosed by elastic membranes in shear flow have shown 
oscillations at frequencies twice the tank-treading frequencies. 

The transition to tumbling behavior at low shear rates occurs when the fluid shear stress 
acting on the particle is no longer sufficient to force the membrane to tank-tread up the 
elastic energy gradient. The surface velocity decreases and the capsule 'solidifies' and begins 
to tumble. In contrast to the tank-trading-to-tumbling transition for a purely viscous drop 
described by Keller and Skalak 2], in which there is no dependence on the shear rate, the 
transition described here depends on the shear rate since f/e ~ 1/7- The transition point is 
a function of the system parameters: f/ec(A, ^). Figure 2 illustrates the tank-treading- 
to-tumbling transition for the geometric (ai/a2 = 4, ai/a^ = 1) and elastic {Eq = 10~^''J) 
parameters of a red blood cell immersed in a viscous fiuid (A = 1.5) being forced at a variable 
shear rate. 

Beyond the critical shear rate the system exhibits type I intermittent behavior 
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i.e. a series of oscillations followed by a tumble followed by another series of oscillations and 
so on. This transition is most easily understood by studying the return map onto the line 
= in the 6-(j) plane modulo vr. For the initial condition 6 = Oq, (p = the trajectory 
crosses the = line again at F{6q), thus defining the return map. Periodic trajectories 
appear as fixed points in the return map F{6q) = 6q and the tank-treading-to-tumbling 
transition occurs via a saddle-node bifurcation in which an unstable limit cycle and a stable 
limit cycle coalesce and disappear. Near the bifurcation Ue — Uec <^ 1 and we expand the 
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return map about 9* ^ where F{9q) — 6q is smallest, so that 

where ^ and i] are constants. Taking 0^ to be the angle after n applications of the map F, we 
see that in the vicinity of the constriction, where F{0o) ^ 0q and ^{0o — 0*)'^ <^ 1, the system 
exhibits nearly periodic oscillations and 0n ~ nrjiUe — Uec)- When ^{0n — 0*Y ~ V^Ue — U^c) 
the system has left the constriction and intermittent behavior occurs. Hence, the number 
of iterations of the map prior to the intermittent behavior scales as 

l/V(f/e-f/ec), 



n ~ 



as seen in figure 4. As Ue is increased further the system enters a regime where the intermit- 
tent behavior consists of a series of tumbles followed by a tank-treading revolution. Finally 
for Ue > Uec* the system fails to exhibit full tank-treading rotations and simply tumbles. 
In this case, the relevant return map is onto the line ^ = in the — (p plane, so that for 
the initial condition = (po, = the trajectory crosses the = line again at -F(0o)- 
This transition is also a type I intermittent transition caused by the coalescence of two limit 
cycles. The phase diagram illustrating the different regimes is shown for both the capsule 
and red blood cell geometries in figure 3. The insets in figure 3 show the return map in 
different regions and two examples of intermittent trajectories. 

By modifying ^ to include visco-elastic membrane properties, our framework can also 
account for the shear flow induced oscillations observed for oil drops coated with a visco- 



elastic protein layer 
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15l |. The details of the behavior discussed in this letter, particularly 
in and around the intermittent regime, are sensitive to parameter values. Therefore, me- 
chanical characteristics of cells and capsules could potentially be deduced from observations 
of such behavior. 
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FIG. 1: (a) Schematic diagram of a tank-treading red blood cell (left) and capsule (right) immersed 
in a fluid of viscosity /x that is being sheared at a rate 7. The viscosity of the fluid inside the cell or 
capsule is /x' so that the viscosity ratio is A = 9 defines the orientation to the shear flow, Vu is 

the surface velocity due to tank-treading, and 4> defines the membrane orientation to the axes. The 
red blood cell and capsule are modeled as ellipsoids of axes Oj. Ue is the dimensionless ratio of the 
resistance of the elastic membrane to deformation to the surface stress applied by the external fluid. 
6 is the phase lag between 9 and the elastic energy, (b) Capsule parameters: axjai = 01/03 = 1.1; 
A = 0; 7 = 5; [/e = 0.21. (c) Red blood cell parameters: 01/02 = 4, ai/03 = 1; A = 1.5; 7 = 5; 

Ue = 0.18. 
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FIG. 2: Red blood cell dynamics. For a red blood cell geometry (ai/a2 = 4, ai/os = 1) and 
elasticity (see text) with a viscosity ratio A = 1.5 a novel transition from tank-treading to 
tumbling is predicted as the shear rate is decreased. For 7 > 'jg^ the membrane tank treads and 
the cell orientation 9 oscillates with amplitude and average that depend on 7. For 7 < j^c the cell 
tumbles and the orientation of the membrane relative to its unperturbed position oscillates. For 
i(f)c < 7 < jdc the cell exhibits intermittent behavior, with a series of tank-treading revolutions 
followed by a tumble or a series of tumbles followed by a revolution. Insets show the dynamical 
regimes for 7 = 2, 1.35, 0.5. 
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FIG. 3: Phase diagram for the red blood ceh (ai/a2 = 4, 01/03 = 1) and capsule (01/02 = 
O'l/o-'i = 1-1) geometries. For a given geometry, the system behavior is determined by Ue, the 
dimensionless ratio of the elastic resistance to tank-treading divided by surface stress from the 
external flow, and A, the ratio of internal to external viscosity. For the red blood cell geometry, 
the insets (A = 1.5, C/g = 0.6557, 0.65) show how the return map {F{9q) gives the value of 9 when 
the trajectory modulo vr beginning at 6* = ^o, = again crosses = 0) changes as the system 
crosses the boundary between tank-treading and intermittent regions. The stable and unstable 
fixed points of the return map coalesce in a saddle-node bifurcation and intermittent behavior 
ensues. Similarly, for the capsule geometry a return map is shown for (p for A = 20, Ue = 0.6, 0.37. 
Both transitions occur for the two geometries investigated. 
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FIG. 4: Type I intermittency: n is the number of tank-treading oscillations per tumble; Ue — Uec 
is the difference between C/g and its value at the boundary between the tank-treading and the 
intermittent regions. The plot is shown for A = 1.5 and a red blood cell geomet ry ( 01/02 = 4, 
o-i/o-s = !)• The -1/2 slope on a log-log plot is indicative of Type 1 intermittency [ij] (see text), 
indicating the coalescence of two limit cycles, which appear as fixed points on the return map of 
figure 3. 
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